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We develop a multi- valued logic for quantum computing for 
use in multi-level quantum systems, and discuss the practical 
advantages of this approach for scaling up a quantum com- 
puter. Generalizing the methods of binary quantum logic, we 
establish that arbitrary unitary operations on any number of 
d-level systems (d > 2) can be decomposed into logic gates 
that operate on only two systems at a time. We show that 
such multi-valued logic gates are experimentally feasible in 
the context of the linear ion trap scheme for quantum com- 
puting. By using d levels in each ion in this scheme, we reduce 
the number of ions needed for a computation by a factor of 
log 2 d. 
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I. INTRODUCTION 

Binary logic gates and Boolean algebra play an impor- 
tant role in classical and quantum theories of computa- 
tion. The unit of memory for binary quantum computa- 
tion is the qu-bit, a quantum system existing in a linear 
superposition of two basis states, labeled |0) and |1). Any 
computation, however large, can be performed using uni- 
versal logic gates that operate on a small, fixed number 
of bits or qubits. fn the quantum case, a unitary trans- 
formation of any number of entangled qubits can be con- 
structed from logic gates that operate on only two qubits 
at a time a result that has no analog in classical 

reversible logic where three-bit gates are needed to simu- 
late all reversible Boolean functions || . We consider the 
extension of universal quantum logic to the multi- valued 
domain, where the unit of memory is the qu-dit |l2) , a 
d-dimensional quantum system with the basis states, |0), 
|1), . . . , \d — 1)- This offers greater flexibility in the stor- 
age and processing of quantum information, and more 
importantly, provides an alternate route to the scaling 
up of quantum computation. 

As in the binary case, a tensor product of many such 
qudits is essential for the efficient storage of information, 
since the number of dimensions in the Hilbert space scales 
exponentially with the number of qudits in the system. 
Allowing d to be arbitrary enables a trade-off between 
the number of qudits making up the quantum computer 
and the number of levels in each qudit. For example, the 
linear ion trap quantum computer || uses only two levels 
in each ion for computing, although additional levels can 
be accessed, and are typically needed, for processing and 
reading out the state of the ion (t]]. By using d compu- 
tational levels in each ion, we reduce the number of ions 



needed for a computation in this scheme by a factor of 
log 2 d, since the Hilbert space of n qudits has the same 
dimensionality as n(log 2 d) qubits, namely d n — 2™ lo S2 d . 
Given the difficulty of trapping and coherently manipu- 
lating a large number of ions in their vibrational ground 
state in this scheme, a reduction in the number of ions 
offered by a multi-valued memory is an advantage. 

A tensor product of qudits is also essential for the ef- 
ficient processing of quantum information. As in the 
binary case, we build unitary transforms on the whole 
system from logic gates that operate within and between 
qudits, creating entangled superpositions, rather than by 
transforming subsets of a non-entangled, unary Hilbert 
space. These elementary multi- valued gates are necessar- 
ily more complex than their binary counterparts, involv- 
ing a controlled transformation of all the levels of each 
qudit. However, a logical network of these gates becomes 
simpler at this expense, invoking a trade-off between the 
complexity of each gate and the number of gates needed 
for a computation [g| . We implement a multi- level gate in 
the linear ion trap scheme by using multiple lasers to ad- 
dress the different transitions in each ion simultaneously. 
In this approach, each multi-level gate takes less time to 
implement than the equivalent binary gate sequence on 
two-level systems, enabling larger computations within 
the decoherence time. 

Quantum computing in multi-level systems is ideally 
described using a multi-valued basis for logic. The in- 
formation stored in a d-level quantum system is funda- 
mentally non-binary in character, since a measurement 
collapses the system to one of these d levels, specifying 
a single value for the qudit, rather than the log 2 d val- 
ues characteristic of a binary representation of the same 
Hilbert space. Moreover, as the d levels in a single qudit 
need not contain any entanglement, two-bit conditional 
logic among these levels is not well-defined, and cannot 
simulate multi-level unitary transforms in practice. By 
contrast, the entanglement between two different d-level 
systems enables conditional two-qudit logic gates, which 
we show to be the elementary operations of multi-valued 
quantum computing. 

Quantum error-correction codes have recently been ex- 
tended to the multi-valued domain, for correcting er- 
rors in a single qudit and multiple qudits 
Fault-tolerant procedures for implementing two-qudit 
and three-qudit analogs of universal binary gates have 
also been developed jfj,|ll|. A proposal for using a cor- 
related photon pair to represent the ternary analog of a 
qubit has been investigated fH} ], but no general scheme 
for implementing multi-valued quantum logic has been 
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proposed. In the following two sections, we derive a set 
of one- and two-qudit gates that are sufficient for univer- 
sal multi-valued computing, and show that these can be 
implemented using multi-level ions in the linear ion trap 
model. 



II. MULTI- VALUED LOGIC GATES 

We review the gates that are universal for binary quan- 
tum logic in a way that facilitates their multi- valued gen- 
eralization. The universal binary gates belong to a family 
of unitary transforms described by three parameters []3| . 
This derives from the fact that up to an overall phase fac- 
tor, any two-dimensional unitary matrix can be written 
as 1 ' 



cos A — e tv sin A 
e^-^sinA e^cosA 



(1) 



expressed in the basis states of a qubit, |0) and |1). The 
parameters A, v, and <f> are usually taken to be irrational 
multiples of 7r and of each other || , since this allows even 
a single gate in Eq. (|l|) to generate all single-qubit trans- 
forms asymptotically by repeated application. However, 
we find it more useful to consider these three parame- 
ters as arbitrary variables in a simulation, with Y 2 rep- 
resenting a family of gates that can be implemented by 
appropriate choice of three physical controls. One of the 
properties of Y% is that it can transform any known state 
of a qubit to |1). Such a transformation, labeled Z 2 , de- 
pends on the coefficients of the state being transformed, 



Z 2 (c ,c 1 ) = F 2 (cos 1 |ci |, arg[c c*] 
co|0)+ Cl |l) ^ |1>, 



arg[cfl) 



(2) 



where |co| + \c\\ = 1. Y 2 also contains the phase gate, 
X 2 , that advances the phase of |1) without affecting |0), 



X 2 (0)=Y 2 (0,0,0) 



|1> 
10) 



e**|l> 
|0>. 



(3) 



Using these two transformation properties of Y%, we can 
show that the two-qubit gates that are universal for quan- 
tum logic take the form, 



r 2 [Y 2 ] = 
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Y 2 



(4) 



acting in the four-dimensional basis of the two qubits. 
The two-dimensional identity 1 2 acts in the basis of |0, 0) 
and 1 0, 1) , and Y 2 acts in the basis of |1,0) and |1,1). 
Taken together, this transforms the second qubit by Y 2 
conditional on the first qubit being in |1). The family of 
gates, r 2 [Y 2 ] , is universal for binary quantum logic in the 



sense that a unitary transform on any number of qubits 
can be simulated by repeated application of these gates 



on no more than two c 
We generalize Eqs 



: ubits at a time. 

(|-||) to the multi- valued case. We 
define Zd as a family of <i-dimensional transforms that 
maps a known single-qudit state to \d — 1), 



Zd(co, ci, . . . , Cd-i) : 
(k)|0) + ci|l) + --- + cd_i|d-l) 1 * \d-l), 



(5) 



where the d complex coefficients, cq, . . . , Cd-i, are nor- 
malized to unity, yielding 2d — 1 real quantities that 
parametrize Zd- As in the binary case, Eq. (Q) does not 
determine Zd uniquely, since it gives the transformation 
of only one of the d states in the basis, namely the super- 
position state with the coefficients, Co, . . . , Cd-i- Since 
it reduces this superposition to a single specified state, 
\d — 1), we may regard Zd as an instance of the quantum 
search algorithm [|l5| , and relate it asymptotically to the 
Walsh-Hadamard and phase transforms used in this algo- 
rithm. Generalizing Eq. (||), we define the d-dimensional 
phase gate Xd as a function of a single parameter, 



Xd(4>) 



\d-l) 1 ► e^\d- 1); 
\p) i-> \p) for p^d-l, 



(6) 



which advances the phase of \d — 1) by <f> without affect- 
ing any other state in the qudit. It turns out that the 
gates, Zd and Xd, are sufficient to simulate all single- 
qudit unitary transforms. We can implement these gates 
by controlling only 2d real parameters, rather than the 
d 2 — 1 that correspond to generalizing Eq. ([[]), thus 
greatly simplifying the physical realization of single-qudit 
gates. If Yd represents either Zd or Xd, then the multi- 
valued analog of T 2 becomes 
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r 2 [F d ] = 











Yi 



(7) 



acting in the cP-dimensional basis of two qudits. The 
identity ld 2 ~d acts on the states, |0, 0), . . . , \d — 2, d — 1), 
and Yd acts on the remaining d states, \d— 1,0), 
\d— l,d— 1). This transforms the second qudit by Yd 
conditional on the first qudit being in \d— 1). We now 
show that such gates are sufficient for constructing arbi- 
trary unitary transforms on any number of qudits. 

Consider an TV-dimensional unitary transform U act- 
ing on n = \og d N qudits. Each state in the computa- 
tional Hilbert space can be written as a tensor product 
of these n qudits, 



fc = 0, 1, 



|A) = 1*1)1*2). . 

.,N-1; fci = 0,l, 



(8) 



, d — 1 for all i, 



2 



where k\ki . . .k n is the base-d representation of k, with 
|fci) denoting the state of the i th qudit. We will use the 
abbreviation, |fci, k%, . . . , k n ), for |fci)|/c2)- ■ • \k n ). Let the 
eigenstates of U be \^f m ) , for m — 1, 2, . . . , JV, with cor- 
responding eigenvalues e l * m . Each such eigenstate can 
be expanded in the computational basis, 

|* TO ) = co|0) + --- + Cjv _i|iV-l) (9) 
= c |0, . . . , 0) + • • • + cjv-i M - 1, • • • ,d - 1), 

where the coefficients are determined by U. Following an 
argument given by Deutsch |l6[ | , we write U as a product 
of JV unitary transforms, each TV-dimensional, that has 
the same eigenstates and eigenvalues as that of U, 

JV 

U = J2 ^ m \*m)(*m\ = WiWa . WN, (10) 

rn—1 



|* m > >-» e**»|* m >; 

I I'm') ^ |*m') for m! ± m. 



(11) 



The problem then reduces to simulating W m for an ar- 
bitrary m. We decompose W m into two transforms that 
are easier to simulate using elementary gates, 



W_7t Y 7 7 — 1 v 7 

m — n 1 7 * 11 1 ' J 1 1 1 1 1 1 1 1 1 ' J ] n ■ 



(12) 



where Z m and X m are the JV-dimensional analogs of Zd 
and Xd- We require only that Z m transform the m th 
eigenstate to — 1), 



Z m (co,ci,..., CAr _i) : |* m ) m- jJV — 1), 



(13) 



which does not determine the transform uniquely, as in 
the case of Zd- We define X m as the transform that 
advances the phase of |JV — 1) by the m th eigenphase, 
leaving all other computational states unchanged, 

, . ( \N-1) ^ e^HTV-l), , . 

|_ \m ) 1— > |m ) for m ^ JV — 1. 

We need to show that W m = Z^ 1 X m Z m satisfies 
Eq. @. First note that 

Z m 1 X m Z m |* m ) = Z m V*™|JV- 1) = e^|* m ). (15) 

For m' 7^ m, the state Z m |W m /) has no projection along 
|AT- 1), 

(JV - l|Z m |tf m ,) = (* ro |Zt,Z m |* ro ,) = (* m |* ro /) = 0, 
which implies that X m has no effect on Z m \^! m r). Hence, 
Zl n X m Z m |* m ,) = Zl n Z in |* m ,) = (16) 



Combining Eqs. (|10j-g_2j) , we see that Z m and X m are 
sufficient to simulate U. For N — d, this implies that Zd 
and Xd contain all single-qudit unitary transforms. In 
the multi-qudit case, we show that Z m and X m can be 



built from the elementary two-qudit gates, ^[Zd] and 
^[-VJ, and their one-qudit counterparts, Zd and Xd- 
We first show this for the n-qudit analog of T 2 , 

T n [Y d ] = Apply Y d to the n th qudit if and only (17) 
if the first n — 1 qudits are in \d — 1), 

where Yd — Zd or Xd- Eq. ( p"7| ) has a matrix repre- 
sentation analogous to Eq. (^), with l d 2 _ d replaced by 
Id^-d- It is easy to see that X m = r n [JV ( j(\I , m )], since 
X m affects only the last computational state, |JV — 1) = 
\d — 1, . . . , d — 1). It is less apparent that Z m is also 
contained in Eq. (p"7|). Z m and Zd are similar in their 
transformation properties in that both take a superpo- 
sition to a single state. However, Zd acts within the 
state space of a single qudit, while Z m transforms the 
Hilbert space of all n qudits. This suggests that Z m 
can be achieved by using T n [Zrf] to target the last qu- 
dit repeatedly, while successively permuting these states 
with the rest of the states in the computational basis. 
First, applying T n [Z d (c N - d , Cjv-i)] to hp m ) reduces 
the superposition of the last d states in Eq. (g) to |JV — 1), 
shown symbolically as 

{\d-l,...,d-l,0),...,\d-l,...,d-l,d-l)} 

= {\N-d),...,\N-l)} |JV-1). (18) 

The superposition of the next d — 1 states in \*& m ) is 
reduced to |JV — 1) by permuting these with the last d 
states but |JV- 1), 



{|JV-2d+l),...,|JV-d-l)} 
«*• {\N-d),...,\N-2)}, 



(19) 



and using Eq. ( |18| ) again. Continuing in this manner, 
successive blocks of d — 1 states in |\& m ) are permuted 
with the last d states but | JV — 1) , and reduced to |JV — 1) 
using r n [Z(j], until the entire JV-dimensional state, |\& m ), 
has been so reduced, completing the simulation of Z m . 
The permutation of states in Eq. ( |l9| ) can also be done 
using r n [Zc;] and r n [Jfd]. To see this, note that a single- 
qudit permutation is already contained in Zd and Xd- In 
particular, if Pd{p, q) denotes the permutation of \p) and 
\q) for p, q = 0, 1, . . . , d — 1, then 

Pd(p,q) = Z d {c , ■ . -,c d -i) X d (ir) Z d (c , . . . ,c d _i), (20) 



On Co 



72 ' Cr ^P'i 



= 0. 



Permuting two states in the n-qudit computational basis, 

\ji,32,---,jn) \k 1 ,k 2 ,...,k n ); (21) 
ji, ki = 0, 1, . . . , d — 1 for all i, 

can be done one qudit at a time, starting with the first 
qudit. The last n — 1 qudits in \ji,j2, ■ ■ ■ ,j n ) are first 
converted to \d— 1) by applying a single-qudit permuta- 
tion, Pd(ji,d— 1), to each qudit i. Then, conditional 



3 



d-1 
d 

















































V 


..J 




LL... d- 


2 _ 








d 


Lil E 


3 d-1 X r 

M_r 


1... 






IU 121 








M d 


2 






LU d-1 





FIG. 1. Construction of r n [Y d ] from r 2 [Y d ], for d > 2. 

on all but the first qudit being in \d — 1), an analog of 
r n [-Pd(ji, fci)] is applied to permute the first qudit from 
\jx) to |fci). The remaining qudits are then restored to 
their original states by the same single-qudit permuta- 
tions, Pd(ji,d— 1), used earlier. This procedure, 

irUPdUud-l): 

\ji,j2,---,jn) >->■ \ji,d- 1, . . . ,d- 1); 

T a [Pd(ji,h)] : 

1> i * d — 1, . . . , d — 1>; 

|fcl, d — 1, . . . ,d- 1) I ^ |fci, J 2 , ■ ■ • , jn), 

is repeated for each of the n qudits, permuting 
\jx, j2, . . . ,j n ) and \k%, &2, ■ ■ • , without affecting any 
other computational state. 

Thus, any n-qudit unitary operator U can be written 
in terms of the logic gates, r n [Frf], for Yd = Zd or Xd- 
We now show that r n [Y"d] can be built from the two- 
qudit gates, r 2 [Yd], of Eq. ([?]). One way of doing this 
is illustrated in Fig. [I] for d > 2. The horizontal lines 
denote the qudits, with solid lines denoting the n com- 
putational qudits and dashed lines denoting additional 
auxiliary qudits that have been initialized to |0). This 
simulation uses r = [~(n — 2)/(d — 2)] auxiliary qudits 
(|~x] means the smallest integer greater than x), where 
(n — 2)/{d — 2) has been assumed for simplicity to be 
an integer in the figure. The vertical lines represent the 
two-qudit conditional gates, originating from the control 
qudit (which is required to be in \d — 1) for the gate to 
apply) and terminating in a box on the target qudit. The 
boxes with two rows p and q represent T%[Pd(pi <?)], the 
conditional permutation of \p) and \q). The box contain- 
ing Yd represents r 2 [Yd], for Yd = Zd or Xj. We want 



the combination of all these gates to implement r„[Yd], 
applying Yd to qudit n if and only if the first n — 1 qudits 
are in \d— 1). 

Reading from left to right in the figure, the first permu- 
tation, T2[Pd(0, 1)], increments auxiliary qudit n+1 from 
|0) to |1) if and only if qudit 1 is in \d— 1). The second 
permutation, r 2 [-P(2(l, 2)], increments qudit n+1 from |1) 
to 1 2) if and only if qudit 2 is in \d — 1), and so on. Con- 
tinuing this way, we see that qudit n+1 reaches \d— 1) if 
and only if all of the first d— 1 computational qudits are in 
\d — 1). This information is then transferred to the sec- 
ond auxiliary qudit, n+2, by the gate T^P^O, 1)], which 
increments qudit n+2 from |0) to |1) provided qudit n+1 
is in \d — 1). This procedure is carried out sequentially 
through all of the computational states, until finally we 
have the auxiliary qudit n+r reaching the state |d — 1) 
(in the case where (n — 2)/(d — 2) is an integer) if and 
only if all of the first n — 1 computational qudits are in 
\d — 1). Controlled by this last qudit, ^[Yd] then acts on 
qudit n, completing the simulation of r n [Yd]. Although 
not shown in the figure, the two-qudit permutation gates 
T2[Pd(p, q)] are re-applied to the auxiliary qudits at the 
end to disentangle them from the computational basis 
and restore them to |0) for re-use. 

This completes the proof that two-qudit gates of the 
form, T2[Zd] and ^[XJ, together with the one-qudit 
gates, Zd and Xd, are universal for quantum computing. 



III. ION TRAP IMPLEMENTATION 

In this section, we discuss one method of implement- 
ing the gates, r 2 [i?d] and r 2 [Xd], in which each qudit is 
represented by a d-level atom. We use the linear ion trap 
scheme for quantum computing, proposed by Cirac and 
Zoller [||, to model the two-qudit interaction. 

The transform Xd does not affect the populations of 
the d states in the qudit, but only changes the phase of 
\d — 1), relative to the other states. Since only one state 
is affected in the process, Xd is effectively the same as 
its binary counterpart, X2, from a physical standpoint. 
We can implement this transform in the atom by cou- 
pling I d — 1) to an auxiliary state in the atom using a 
27r-pulse, which does not leave any population in this 
state at the end of the pulse. In the interaction picture, 
the phase of \d — 1) after the pulse will be different if the 
detuning is made time-dependent. One way to realize 
a time-dependent detuning is by using a Stark field to 
shift the energies of the two levels over time. The other 
computational states in the atom are not affected in the 
process if they are far off-resonance and the fields are 
sufficiently weak. 

Unlike Xd, the transform Zd involves all of the states 
in the qudit, acting on a (i-state superposition with the 
coefficients, Co, c\, . . . , Cd-x, as shown in Eq. (^|). The 
implementation of such a transform can be posed as a 
problem in quantum optimal control fll7|. If there are 
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d-1 lasers 



FIG. 2. Linear ion trap, with q ions. Trap axis is along x. 



d-1) 




FIG. 3. Level scheme for a d-level ion, with d— 1 neighbor- 
ing transitions. 

2d — 1 physical controls available for manipulating the 
qudit, an optimization on these controls can be done with 
the fidelity governed by Eq. (||). 

A time-domain approach to this problem in a multi- 
level atom was studied by Noel and Stroud [Q , where the 
control parameters were the amplitudes and time delays 
of a sequence of d laser pulses, and the goal was to excite a 
d-state wave packet in the atom, starting from the ground 
state. However, when the ground state population be- 
comes significantly depleted, non-iterative methods may 
not be sufficient for creating arbitrary wave packets in the 
atom ]ig|| . We consider an alternate frequency-domain 
approach to implementing Zd, where the control parame- 
ters are the amplitudes and phases of different laser fields 
that are tuned near resonance to d— 1 atomic transitions, 
and that are adiabatically turned on and off over a con- 
trolled time period. To realize the two-qudit gates T 2 , we 
consider a multi-valued extension of the linear ion trap 
scheme. 

Consider q identical (i-level ions confined in a linear 
harmonic trap with frequency v x , each of which can in- 
teract with d — 1 lasers at a given time (see Fig. |]). If 
each laser is detuned from the associated atomic reso- 
nance by u X) only the center-of-mass (CM) normal mode 
of the trap is excited in the absence of power broaden- 
ing. By laser cooling, the ions are initially assumed to 
be in the vibrational ground state of this mode, where 
each ion vibrates about its equilibrium position with an 
amplitude that is small compared to an optical wave- 
length. The trap is then characterized by a Lamb-Dicke 
parameter, r\ = k x {hj2mv x ) l l' 1 , that is small compared 
to unity, where m is the mass of each ion and k x is the 
laser wave vector along the trap axis. If d^ and a are the 
creation and annihilation operators for the CM mode, 
and ijjj — are the internal projection operators for 



a given d-level ion in the trap, the Hamiltonian for this 
ion in the absence of interaction fields is 

d-1 

Hq = hi/ x (d J< d + 5) + hujj ajj. (22) 
3=0 

The computational level scheme considered is shown 
in Fig. |^, where the transition frequencies, tOjj+i = 
\ujj+i — are distinct compared to the linewidths of 
the levels. For the purpose of implementing Z e i, it is suf- 
ficient to have only the neighboring levels coupled. This 
can be reinforced by using appropriate selection rules to 
suppress the other transitions. The d — 1 neighboring 
transitions are driven by near-resonant laser fields that 
have a standing- wave configuration along the trap axis, 

E(i,t) 

d-2 

= ^2 e j j'+i l E j,j+i e~ M " +lt + c.c] cos(kjj + ix + tp) 

(23) 

d-2 

= X)ew+i[^, j+ ie- fa «+ 14 + c.c] 

J x [ cos(¥;) _!^±l ( at + a) S in^) + o( v l j+1 )}, (24) 

where Ejj+i and otjj+i are the (complex) field ampli- 
tudes and field frequencies corresponding to the atomic 
transitions, and and are the associated po- 

larizations and wave vector components. The field depen- 
dence on y and z has been suppressed due to the strong 
trap confinement along these directions. When ip = tt/2 
or 0, the standing waves make a node or antinode at the 
ion's equilibrium position, (x) = 0. We have used 

x = {h/2qmv x ) 1/2 {a) + a) (25) 

for the displacement of the ion from equilibrium, where 
qm is the effective mass of the CM mode. Each cosine in 
Eq. ( p3| ) has been expanded in powers of the correspond- 
ing Lamb-Dicke parameter, ffoj+i = kjj + i(h/2mi' x ) 1 ^ 2 , 
and only terms up to first order in ?7j,j+i are kept in the 
limit, 77j,j+i -C 1, for each j. For the level scheme under 
consideration, the internal dipole moment of the ion is 
effectively 

d-2 

= [ d «+i4j+i + d jj+i^j+i]' ( 26 ) 

where (Jjj+i and d* J+1 are the transition operator and 
matrix element corresponding to the downward transi- 
tion between levels \j) and \j + 1). The ion-field interac- 
tion is described in the dipole approximation, using the 
Hamiltonian, 

H dip = -d-E{x,t), (27) 
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where the field depends on the center-of-mass position, x, 
of the ion. This interaction couples the electronic (a, <t') 
and vibrational (d,a^) degrees of freedom of the ion. 

We study the time evolution in the interaction picture, 
where the operators a(t) and djj+i(t) evolve according 
to H , as ae^-' and &j,j±i e~ Wj -' +lt . When ip = 0, 
the field expansion in Eq. (E4) does not contain and a 
to first approximation, and i/dip only affects the internal 
states of the ion. Tuning each laser to resonance in this 
case, — uij.j+i, we find that ifdi p becomes time- 

independent in the interaction picture under the rotating- 
wave approximation, 



cZ-2 
3=0 



^'j+i + ^ij+i ^iJ+i]' 



(28) 



where f^j.j+i = (dj,j+i • ej_j + iEjj + i)/h is the Rabi fre- 
quency for the transition between levels \ j) and \j + 1). 
Alternately, when tp = tt/2, the field is linear in and a 
to first approximation, and H^ip affects both the internal 
and external states of the ion. In this case, detuning each 
laser above or below resonance by the trap frequency, 



z 3,3 + 



i = ujj j + i ± z/ x , we find 



d— 2 



-ftdip,U + = ^ 



\ - Vj.j+i 



+ / y 

3=0 



(29) 



H, 



dip,U_ 



(30) 



The unitary time evolution operator corresponding to 
Eq. (§), 



V = exp[-i(t/h) #di P ,v], 



(31) 



mixes the d internal states of the ion without affecting 
the trap state, and turns out to be sufficient for gener- 
ating the single-qudit gates Zd, up to an overall phase 
factor. The time evolution operators corresponding to 
Eqs. §§, 



U± = exp[-i(t/H) Hdip.uJ, 



(32) 



conditionally couple the internal and external co- 
ordinates of the ion. Whenever the internal energy of 
the ion is raised (<5^), J7+ raises the trap energy (d^), 
while i/_ lowers the trap energy (d), in tandem. This 
conditionality arises from the rotating-wave approxima- 
tion, which retains only the energy-conserving terms in 
the Hamiltonian. Using U± and V in stages, we will show 
that the two-qudit gates, r 2 [ld], can be implemented be- 
tween two ions in the trap. 

First we show how to construct Zd from V. In the 
binary case, d = 2, we set all the Rabi frequencies except 



Slo,i equal to zero. The levels |0) and |1) then undergo 
two-level Rabi oscillations, 



V 



ti- 



ne iQ^S 
in 0A s n'c 



|o) 
ID 



(33) 



where C — cos fit, S — sin fit, and tt = |f2o.i|- Given a 
state cq|0) + ci|l), we can choose J7o,i and t such that 



ilcocil 



cos Vtt 



(34) 



which makes V implement the transform of Eq. (^|), up 
to an overall phase, iargci, simulating the binary gate 
Z 2 . In the ternary case, we set all the Rabi frequencies 
except f2o,i and fiio equal to zero, leaving a three-level 
A-system (see Fig. ||). In this case, the levels |0), |I) and 
1 2) evolve according to 



V 



|^i,2| 2 + |fio,i| 2 c ifi* 01 ns ^iOi, 2 (c-i) 



(35) 
|0) 
|2) 
\^ 2 \ 2 - 

Given a state co|0) + ci|l) + c 2 |2), we can choose fio,i> 
Oi 2 and t such that 



n 0:1 fi* 12 (c - 1) ifi* 12 ns \n 0A \ 2 + \n h2 \ 2 c 



where C = cosilt, S — sin fit, and £1 2 = |f2o,i| 



0,1 



V. 



ilcil 2 



1-C 



n 



(36) 



cos ilt = 



1- C 2 



- 1, 



which makes V implement the transform of Eq. (|^) for 
d = 3, up to an overall phase, iargC2, simulating the 
ternary gate Zj,. In the d- valued case, we require V to 
implement in Eq. (|^) for an arbitrary d, 



1^(^0,1, fil, 2, ■ ■ • , ^d-2,d-i; i) : 
co|0)+ci|l) + --- + c d _i|d-l) i ^ e**|d- 



1) 



(37) 



where <fi has been introduced to allow for an overall phase 
offset. The controls in V are the d — 1 complex Rabi 
frequencies, and the interaction time. The adjoint of the 
transform in Eq. (|37j ) can be written as 

V*\d-1) = e-** [co|0) + ci|l) + --- + c_i|d-l)]. 

(38) 

Projecting this equation onto (p\ and taking the complex 
conjugate of both sides, we get 



(d-l\V\p) 



c v e , P 



0,1, 



1. 



(39) 



which relates d of the matrix elements of V to the co- 
efficients, Co, . . . , Cd-i- These d equations have to be 
inverted to find the controls, fio,i> ■ ■ ■ > ^d— 2,rf— lj and t. 
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Analytical solutions are given for the binary and ternary 
cases in Eqs. (|34|) and ([56]). This allows V to implement 
Zd up to a phase gate, Xd{4>)- 

Two-qudit gates of the form ^[Id] can be imple- 
mented using both U± and V interactions, and an aux- 
iliary manifold of d additional levels in each ion. To see 
this, write the original state of the two-ion system in the 
form, 



I«-)cI $ )tIQ), 



(40) 



where |\&) c is the original control ion state, |3>) T is the 
original target ion state, and |0) is the trap ground state. 
Applying a 7r-pulse of the U± interaction to the control 
ion, we first transform all of the computational states 
in this ion except \d— l) c to their auxiliary counter- 
parts, conditional on exciting the trap to We leave 
\d— l) c |0) unaffected by turning off the corresponding- 
laser in U±. We then restore the internal state of the 
control ion to its original configuration by using a 7r-pulse 
of the V interaction, which does not affect the trap. The 
entangled state of the system is then given by 



|*d-l> C l*> T IS) + l*other)cl*> T ID. 



(41) 



where |\P) C = \^d-i)c + l*othcr) c - We see that all of the 
control states except \d — l) c are entangled with the trap 
state |1). Applying a 7r-pulse of the U± interaction to 
the target ion now, we transform all of its computational 
states to their auxiliary counterparts, conditional on de- 
exciting the trap, 



l*d-i; 



1$) 



|*otho 



I*. 



(42) 



where |$ a ux)T is the original target ion state written in 
the auxiliary basis. The first term in expression (^l|) is 
not affected by this operation since the trap ground state 
cannot be de-excited. Next, applying V in the compu- 
tational basis of the target ion, we simulate Yd = Zd or 
Xd, transforming |$) T but not affecting |$ a 



|*«,_i> c {Y d \$) T } |0) 



$aux) T |Q). (43) 

is then restored to the com- 



| Mother) C 



The target ion state l^aux^ 
putational basis by reversing the operations that took us 
from expression (HG) to expression (E2J) , giving 



\^d-i) c {Y d \^} T }\0} 



|*othcr> C |$> T IO>. 



(44) 



This completes the implementation of T2 [Yd] on expres- 
sion ([40|), with the target ion transformed by Y d condi- 
tional on the control ion being in \d— l) c . This two- 
qudit logic is made possible by the U± interaction, which 
allows the information about whether the control ion is 
in \tyd-i)Q or | Mother) c to be carried to the target ion via 
entanglement with the trap states, |0) or |1). This shows 
that universal multi-valued computing is feasible in the 
linear ion trap scheme. 



IV. SUMMARY 

We conclude with a comparison of the binary and 
multi-valued approaches to quantum computing. The 
main advantage of the latter is a logarithmic reduction in 
the number of separate quantum systems needed to span 
the quantum memory. For a Hilbert space of N dimen- 
sions, corresponding to n 2 = log 2 N qubits, the number 
of qudits needed to store this information is 



log 2 ^ 
log 2 d 



n 2 
\og 2 d 



(45) 



Note that this retains the same scaling in N and d with 
the inclusion of the auxiliary qudits used in the gate con- 
struction of Fig. [I]. Using a binary equivalent of this 
construction, we find that the overall time-complexity of 
a binary simulation is 0[n 2 N 2 ]. That is, this many two- 
qubit gates are required to simulate an A-dimensional 
unitary operator U. By analogy, the number of two- 
qudit gates used in the construction of section || scales 
as 0[n 2 N 2 }, or 



0[n 2 N 2 ] = O 



(log 2 N) 2 N 2 



(log 2 rf) 2 



O 



n 2 N 2 



(log 2 d) 5 



(46) 



where we have used Eq. (45). Eq. ( |46| ) represents an up- 
per bound on the time-complexity of a multi- valued sim- 
ulation, and shows that this has a (log 2 d) 2 advantage 
over the binary case. This comes at the cost of larger 
elementary gates, ^[Yd], which require d states in each 
ion to be controlled, not just two. This suggests that 
we ought to multiply the multi-valued time-complexity 
in Eq. ([46| ) by d for a physically relevant comparison 
with the binary case. However, the frequency-domain ap- 
proach taken to constructing the V and U± interactions 
in Eqs. (p8[-30|) assumes that the d — 1 lasers operate si- 
multaneously on each ion, which allows a two-qudit gate 
to be implemented without slowing down the computa- 
tion in real time. The cost of the multi- valued speed-up 
in this case is the need for multiple lasers to address the 
corresponding transitions in the d-level ion. 

Finally, we must make note of the non-logarithmic scal- 
ing in N in Eq. (|4^), which shows the inefficiency of this 
construction for simulating arbitrary TV-dimensional uni- 
tary transforms. This is analogous to the binary case, 
where only certain unitary transforms admit an efficient 
simulation in terms of elementary gates, making them 
useful for efficient quantum algorithms. 

We thank David Aronstein for helpful comments. 
This work was supported by the Army Research Office 
through the MURI Center for Quantum Information. 
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